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We describe how to overcome some problems that usually prevent from obtaining an efficient algorithm to fix 
a generic covariant gauge on the lattice. This gauge is the lattice equivalent of the generic gauge usually adopted 
in perturbative calculations. It depends on a parameter whose value can be varied in order to check the gauge 
dependence of measured matrix elements. 
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We discuss how to enforce on the lattice the 
generic covariant gauge: 



9^A^«(")(a;)=A(a;); 



fl) 



where A{x) are matrices belonging to the Lie al- 
gebra of the SU(3) group and a is the gauge pa- 
rameter which is the width of the gaussian distri- 
bution used to extract the A(a;) matrices as will 
be described in the following. The differential 
equation (|I|) is supplied with periodic boundary 
conditions on the lattice. 

Before going into details, we want to make a 
preliminary remark concerning the definition of 
the gauge potential Af^ on the lattice. As it is 
well known the gauge potential, on the lattice, 
is a derived quantity being the links C/^ the nat- 
ural fields. The link is the parallel transporter 
of the theory and is connected with the gauge 
potential A^^ only in the continuum limit. Of 
course, it is necessary to choose a definition of 
Afj_ in terms of the links in order to fix the gauge; 
this is not a problem because any A^ definition 
differing from another by irrelevant terms ( 0{a)) 
must go to the same continuum operator up to a 
constant. This trivial, though fundamental, con- 
sideration was checked non-perturbatively by nu- 
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merical simulations, H] obtaining a spectacular 
coincidence (see Fig. ^ between the rescaled val- 
ues of the correlator: 



3^2 



I x,y 



built up assuming the following different A^ def- 
initions: 



M^) 



{U^ix)-UUx)) 
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Af,{x) = 
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{(U^JX))^ - {Ul{x)f)traceless 
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We found that the expected relation A' Jx) = 
C{go)A^,{x) + O(a^) is fulfilled and the effect of 
order a terms is negligible in the measured ma- 
trix elements. This result substantially supports 
the correctness of the continuum limit of matrix 
elements built up by the insertion of operators 
representing ^^ on the lattice. 

Actually the freedom of the A^ definition on 
the lattice, instead of being considered problem- 
atic, has been rarely exploited in literature and 
it is interesting to quote two recent papers that 
take into account this possibility for different pur- 
poses. Ref. [g proposed to use, in the Landau 
gauge, the exponential form for the relation be- 
tween C/p and Af^ instead of the linear one (^. 




Figure 1. Comparison of tlie matrix elements of 
(AiAi){t) (crosses) and tiie rescaled {AiAi) ■ Cf{go) 
(open circles) as function of time for a set of 50 ther- 
malized SU{3) configurations at /3 = 6.0 with a vol- 
ume 1/ ■ r = 8^ • 16, Ci = .729(2). The data have 
been slightly displaced in t to help eye, the errors are 
jacknife. This figure is taken from ref. |lf. 



In ref. B it has been attempted to "improve" 
the Landau gauge fixing by studying different A^ 
definitions. 

The classical procedure to fix a generic covari- 
ant gauge on the lattice can be read from the 
following expression of the expectation value of a 
gauge dependent operator O on the lattice: 
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G{a))^-J^ J d-'xTriA^J^-^fSSiU) 
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where S{U) is the Wilson lattice gauge invariant 
action, and G{a) is the gauge transformation that 
enforces the gauge condition, eq. (|l|). 

Following the usual procedure, the gauge trans- 
formation G{x) must be determined numerically 



by applying a steepest descent iterative algorithm 
which minimizes an ad hoc functional F[G] cho- 
sen in order to have minimum (extremum) cor- 
responding to a gauge transformation satisfying 
the gauge condition. 

The definition of < O > shows that before 
applying to each thermalized configuration the 
gauge fixing algorithm one has to generate a set 
of A, one for each lattice site, distributed as a 
gaussian of width a. Instead of taking more dif- 
ferent sets of A as it is suggested by (||) and then 
averaging, it is reasonable to take only one set of 
A for each configuration. In fact one generates for 
each thermalized configuration a different sets of 
A matrices with the same width a. 

The functional -F[G], that in the Lan- 
dau gauge is given by Fi[G] = ||A'^|P = 
/Tr (^A^{x)A'j^{x)) d'^x, must be changed in or- 
der to reach the new gauge condition. This prob- 
lem was discussed in B where the following more 
complicated functional form was suggested: 
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(9^A^(") - A)^ 
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The stationary points of this functional corre- 
spond to the following gauge condition: 



C.5,(5,,A;;(")-A) = 



(7) 



Of course, when S^Ap 



G(a) 



(x) — A{x) = 0, the new 
functional is stationary but there can be other 
non trivial solutions. This fact is obviously a 
problem of the new functional, and simpler func- 
tional forms are difficult to find, as it has been 
discussed in Q. We will assume the functional 
form (pi) and we will ignore the problem of spuri- 
ous solutions in the following. The choice of the 
form of F[G] is not sufficient to solve the numer- 
ical problem of fixing the generic covariant gauge 
because, due to the comphcated structure of (0), 
it is likely that the usual gauge fixing algorithm 
either does not converge or takes too much com- 
puter time. This problematic convergence can be 
traced back to the following reason. In the Lan- 
dau case, being Fl proportional to U^, when the 
gauge fixing algorithm visits the lattice point x it 
puts only the gauge transformation in the point 
X different to the identity and then locally Fl de- 
pends linearly on G{x). Of course this property 



is not valid anymore for a naive transcription on 
the lattice of F[G] because this shows a natural 
quadratic dependence on G{x). 

At this point one can take advantage of the 
freedom to choose the discretisation of A^. In 
fact it is possible to find a discretisation ( "driven 
discretisation") for the functional F[G] such that 
it takes only a local linear dependence on G{x). 
This aim can be reached by allowing different 
terms of F to be chosen in such a way to guaran- 
tee the local linear dependence on G{x), instead 
of being an algebraic consequence of a particular 
definition of A^. In the following we show some 
examples of the discretisation for the terms con- 
tained in the functional F[G]: 

2iagA^{x) = U^{x) - Ul{x); 
a^g^Af,{x)A^{x) = 

I + Ul{x)U,{x)-UAx)-Ul{x)- 
a^g^A^{x — ii)Ai,{x — v) = 1+ 

Uf_,{x - ti)Ul{x -v)- Uf_,{x - n) - Ul{x - v). 

Each term locally transforms linearly in G{x) and 
it is a possible discretisation of the continuum op- 
erator written on the left side of each expression. 
Using the driven discretisation the expression of 
F[G\ assumes the following compact form: 



where 

J{x) = 



-81 



[Ul{x -v) + Ui,{x)] - iagA 



As it is shown in Fig. pi the number of sweeps 
necessary in order to reach a prefixed quality of 
the gauge fixing increases when a decreases. Each 
sweep takes around 20% more computer time 
than in the case of the functional Fl ■ The qual- 
ity of the gauge fixing is measured by adopting as 
quality factor the variation SF of F[G] between a 
sweep and the next. 

Further details on the formalism, spurious so- 
lutions and checks will be discussed in a forth- 
coming paper (Ref. [Q). 




Figure 2. Behaviour of 5F vs gauge fixing 
sweeps (with overrelaxation) for different values of a; 
(3 — 5.8, volume=6'*. 
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